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of  a limited  amount  of  information.  By  considering  the 
mathematical  relationships  between  the  detectors,  formulations 
have  been  developed  in  this  dissertation  which  allow  for  the 
assessment  of  both  the  amount  and  associated  spatial 
uncertainty  of  the  radioactive  material  measured  from  N 
detector  responses. 

The  analysis  of  the  proportional  relationships  between 
the  detector  responses  results  in  the  following:  the 
calculation  of  the  possible  activities  from  a detector 
response  of  unknown  activity  and  distribution,  the  calculation 
of  the  largest  possible  uncertainty  from  N measurements  for  a 
constrained  distribution  of  activity,  the  property  that  local 
optimization  over  an  unconstrained  distribution  is  a global 
optimization.  The  analysis  of  the  first  and  second  derivatives 
of  a response  sets  results  in  a technique  for  reducing  both 
the  number  of  optimization  steps  and  the  number  of  independent 
variables  necessary  for  the  calculation  of  the  largest 
uncertainty. 

Increasing  the  number  of  symmetric  detectors  surrounding 
a container  results  in  the  largest  uncertainty  approaching 
some  asymptotic  limit.  A relationship  between  an  infinite 
number  of  symmetric  measurements  and  a definite  integral  is 
developed  allowing  the  solution  to  the  asymptotic  limit  to  be 
obtained  by  analytical  means. 

The  optimization  of  a detector  system  by  minimizing  the 
largest  uncertainty  with  respect  to  detector  position  is 


Utilizing 


CHAPTER  1 
INTRODUCTION 


The  focus  of  this  dissertation  is  the  application  of 
dimensional  space  for  the  assay  of  randomly  distributed 
radioactive  material.  The  presence  of  radioactive  material  can 
be  determined  by  performing  measurements  with  suitable 
detectors  (Knoll  1979) . The  utilization  of  these  detector 
responses  for  the  assessment  of  a measured  amount  of 
radioactive  material  is  well  understood  for  known  source 
distributions  of  the  radioactive  material.  However, 
uncertainty  concerning  the  radioactive  material  distribution, 
possibly  due  to  some  inherent  random  property  or  lack  of 
sufficient  information  on  the  part  of  the  assessor,  poses 
additional  complications. 

Given  the  assumptions  that  the  spatial  distribution  of 
the  radioactive  material  within  the  bounds  of  a container  is 
unknown  and  that  a limited  number  of  radiation  detectors  are 
present,  it  is  difficult  and  often  impossible  to  construct  the 
actual  spatial  configuration  of  the  source  present.  In  many 
cases,  different  amounts  of  radioactive  material  with 
different  spatial  distributions  can  result  in  identical 
detector  responses.  The  methodology  presented  in  this 


dissertation  is  applicable  to  the  problem  of  effectively 
utilizing  the  limited  information  container  from  a relatively 
small  number  of  detector  measurements  and  not  on  the  more 
sophisticated  problem  of  mapping  the  concentrations  of  the 
radioactivity  present. 

The  introduction  of  mathematical  methods  for  the  assay  of 
randomly  distributed  radioactive  material  has  been  limited. 
The  difficulties  of  assessing  both  the  amount  and  the  inherent 
uncertainties  of  the  radioactive  material  in  the  presence  of 
random  or  unknown  distributions  has  been  suggested  as  an  area 
requiring  prior  information  and  expert  judgment  (Sampson 
1992) . If  the  distribution  of  the  radioactive  material  can  be 
assumed,  probabilistic  methods  can  be  applied  (Ben-Haim  1983) . 
Techniques  such  as  rotation  of  the  source  container, 
segmentation  of  the  measurement  volume  seen  by  the  detector 
and  linear  combinations  of  detector  responses  have  been 
introduced  as  methods  for  reducing  the  assessment  uncertainty 
for  such  configurations  as  a point  source  and  a small  cylinder 
source  within  a waste  container  (Martens  and  Fliss  1990) . 

The  application  of  dimensional  space  has  previously  been 
introduced  as  a tool  for  a calculating  a performance  criterion 
called  the  relative  mass  resolution  (RMR)  . The  RMR  is  defined 
as  the  ratio  of  the  maximum  amount  of  radioactivity  compared 
to  a fixed  amount  of  activity  that  could  result  in  an 
identical  set  of  detector  responses.  This  application, 
however,  does  not  extend  to  the  problems  of  either 


interpreting  multiple  detector  responses  in  the  presence  of 
randomly  distributed  radioactive  material  (Ben-Haim  1983)  or 
assessing  the  actual  uncertainty  from  a given  set  of  detector 
measurements.  In  addition,  the  search  method  suggested  for  the 


imber  of  limitations.  The  search  method  is  a 
local  optimization  method  (Scott  1991)  and  is  valid  only  for 
cases  where  the  detector  response  set  is  convex  (Ben-Haim 
1985a,  1985b) . Appendix  A contains  a number  of  definitions 
concerning  tf-dimensional  space  and  convex  set. 

The  above  limitations  are  addressed  within  this 
dissertation  by  extending  the  calculational  methods.  A 
solution  to  the  problem  of  interpreting  multiple  detector 
measurements  in  the  presence  of  spatial  uncertainty  is 
proposed.  A global  search  technique  for  convex  sets  is 
presented.  A complementary  derivative  technique,  which  reduces 
the  number  of  search  variables,  is  added  to  the  existing 
method  for  calculating  the  RMR  (Ben-Haim  1983) . A method  of 
analyzing  the  asymptotic  limit  of  an  increasing  number  of 
detectors  is  presented. 

Incorporated  within  the  development  of  the  RMR  is  the 
concept  of  representing  the  set  of  detector  responses  from 
unconstrained  source  distributions  by  the  set  of  detector 
responses  from  point  source  distributions.  A number  of 
examples  and  proofs  expanding  on  the  subject  of  convex  sets 
exist  (Ben-Haim  1985b)  . A related  dissertation  (Chan  1988)  was 


;h  relied  heavily 


formulated  concepts  and  equations. 


on  previously 


while  unified  by  the  concept  of  dimensional  space,  has  taken 
a number  of  directions.  Published  literature  in  the  field 
(Ben-Haim  1985b)  mentions  a number  questions  and  problems  that 

research  were  briefly  discussed.  Solutions  to  two  of  the 

problem  (adaptive  assay-advanced  concepts)  are  given  in  the 
Appendix  B.  Along  with  the  unresolved  issues  that  can  be  found 
in  the  published  literature,  the  practioner  of  the  assay  of 
spatially  random  material  is  presented  with  a number  of 


variables  as  the  number  of  measurements  increases, 
investigations  have  been  performed  in  the  area  of  variable 
reduction  and  optimization  of  search  techniques. 

The  structure  of  the  dissertation  is  divided  into 
methodologies.  chapter  3 extends  the  application  of 
dimensional  space  to  the  problem  of  interpreting  detector 
responses.  The  resulting  method,  called  the  ray  method,  can  be 
utilized  to  assess  the  amounts  of  radioactive  material  that 
could  result  in  a given  N detector  response.  Chapter  4 extends 
the  ability  to  calculate  the  RMR.  The  use  of  supporting  planes 
was  previously  developed  (Ben-Haim  1983)  as  a method  for  the 
calculation  of  the  RMR  of  a convex  response  set  utilizing  only 
the  fundamental  set  that  generated  the  convex  set.  The 
contributions  to  this  method  involve  the  ability  to 
efficiently  calculate  the  RMR.  For  some  problems  the  weight 
coefficients  can  be  substituted  for  derivatives,  thus  reducing 
the  number  of  variables  of  the  search  space.  In  addition,  the 
difficulties  of  utilizing  the  existing  search  method  for 
calculating  the  RMR  is  discussed  and  alternative  strategies 
are  developed.  The  ray  method  is  applied  to  the  solution  of 
the  RMR  in  the  presence  of  nonconvex  sets  and  to  the 
development  of  a global  optimization  method  for  the  solution 
of  the  RMR  in  the  presence  on  convex  sets.  Chapter  5 develops 
a method  for  determining  the  asymptotic  limit  of  the  RMR  as 
the  number  of  symmetric  detectors  approach  infinity  in  the 
presence  of  a fixed  number  of  nonsymmetric  detectors. 


fundamental 


theory  that  have  previously  been  developed  (Ben-Haim  1985a) . 


A coordinate  system  is  constructed  with  the  x-axis 
representing  f,,  the  y-axis  t,  and  the  z-axis  fj.  A response 


of  source  activity.  Figure  1-2  represents  the  above  detector 
given  source  distribution. 


Assuming  a constant  attenuation  coefficient,  if  the 
location  of  a point  source  of  radioactivity  continuously 


The  application  of  much  of  the  methodology  presented  is 

proportional  to  the  amount  of  activity  present  and  that  the 
self-attenuation  of  the  source  is  negligible.  To  allow  for 
analytical  solutions  in  the  examples  presented,  the 


simplified  by 


assumptions.  The  detectors  respond  as  point  detectors. 
Attenuation  of  both  the  source  container  and  the  medium 
between  the  container  and  the  detector  is  considered 
negligible.  The  response  from  the  nth  detector  to  a point 
source  is  given  as  where  d,  is  the  distance  from  the 
point  source  to  the  detector. 

While  the  main  purpose  of  the  dissertation  is  to  extend 
the  application  of  dimensional  space  to  the  problems  of 
radioactive  material  assessment,  a secondary  purpose  is  to 
increase  the  understanding  of  the  subject  for  the 
practitioner.  The  utilization  of  dimensional  space  for  the 
assay  of  radioactive  material  is  not  common.  In  addition,  the 
emphasis  of  this  dissertation  is  on  the  understanding  of  the 
mathematical  tools  possible  from  applying  dimensional  space 
rather  than  on  the  utilization  of  supporting  planes  or  convex 
sets.  The  result  of  such  an  approach  is  that  while  a number  of 
innovative  ideas  are  presented  from  a unified  perspective, 
significant  work  concerning  concepts  previously  developed, 
such  as  the  RMR  and  convex  sets  (Ben-Haim  1984,  1985a)  has 
been  summarized. 

A number  of  examples  are  presented  which  illustrate  the 
ideas  and  the  problem  solving  technigues  formulated.  For 
graphical  reasons  two  detectors  are  usually  stipulated. 
Several  detector,  source  container  configurations  are  repeated 
throughout  the  examples.  The  one  dimensional  source  container 
presented  is  a line  segment  on  the  x-axis  from  x=-0.5  to  x=0.5 
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CHAPTER  2 

FUNDAMENTAL  CONCEPTS 

Application  of  Vector  Space  for  Detector  Responses 

the  amount  of  radioactive  material  present  and  that  the 
radioactive  material  distribution  is  known,  the  use  of 
additional  detectors  at  different  spatial  locations  does  not 

to  calculate  the  amount  of  radioactivity  present.  A common 
example  is  the  use  of  a single  detector  response  at  a specific 


thus  allowing  a calculated  value  without  uncertainty  to 
represent  the  activity  present. 

in  uncertainty  concerning  the  interpretation  of  a detector 
response.  When  attempting  to  assess  an  amount  of  radioactivity 
present  as  illustrated  in  Figure  2-1,  one  unit  of  activity  at 
a specified  location  could  result  in  the  same  response  as  two 


demonstrates  the  ability  of  a two  detector  system  to 
distinguish  between  these  two  distinct  cases  of  source 
activity  and  location. 

The  detector  system  in  Figure  2-2  can  be  described  as  a 
two  variable  array  (f[,  fj) , where  f,  is  the  response  of  the 
first  detector  and  f2  is  the  response  of  the  second  detector. 
The  utilization  of  the  two  detector  system  in  Figure  2-2 
results  in  a less  uncertain  assessment  of  activity  than  the 
one  detector  system  illustrated  in  Figure  2-1. 

If  the  spatial  distribution  of  the  activity  is  known,  the 
methodology  for  assessing  the  radioactive  material  for  a one 
detector  system  with  a known  distribution  can  be  performed  by 
utilizing  a previously  determined  response  per  activity  factor 
for  a particular  nuclide.  The  activity  is  calculated  by 
dividing  the  factor  into  the  given  detector  response. 

For  a one  detector  system,  if  the  spatial  distribution  of 
the  activity  is  unknown,  then  the  maximum  and  minimum  possible 
detector  responses  for  one  unit  of  activity  are  calculated. 
For  the  detector  system  in  Figure  2-1,  the  cases  of  greatest 
and  least  detector  response  are  point  sources  located  at  the 
distances  closest  to  and  farthest  from  the  detector.  The 
amount  of  activity  possible  from  a given  detector  response  is 
calculated  by  dividing  the  given  detector  response  by  the 
minimum  and  maximum  detector  responses  to  obtain  the 
respective  maximum  and  minimum  amounts  of  activity.  All  other 
possible  amounts  of  activity  that  could  result  from  the 


extreme  detector  responses  are  between  these  calculated 
amounts  of  activity. 

When  attempting  to  apply  the  above  one  detector 
intuitive  type  procedure  to  a multiple  detector  system 
additional  questions  arise.  What  is  an  extreme  detector 
response?  How  is  a response  "between"  to  other  extreme 
responses  defined?  How  is  an  assessment  of  activity  performed 
from  a multiple  detector  response?  The  above  questions 
indicate  the  failure  of  attempting  to  use  a strictly  intuitive 
method  for  multiple  detector  responses.  A more  mathematical 
approach  however,  is  possible  and  has  been  previously 
presented  (Ben-Haim  1983) . From  this  approach,  which  applies 
N dimensional  space  to  an  N detector  response  system,  the 
concepts  of  extreme  responses  and  responses  between  responses 
are  easily  extended  from  a one  detector  system  to  a multiple 
detector  system. 

The  foundation  of  the  mathematical  approach  is  that 
responses  from  N different  detectors  are  considered  as 
coordinate  axes  in  S dimensional  space.  The  application  of  N 
dimensional  space  allows  for  the  representation  of  detector 
responses  as  points.  A set  of  points  resulting  from  a specific 
amount  of  activity  is  called  a response  set.  An  extreme  point 
is  the  point  of  farthest  distance  from  the  other  points.  The 
property  of  proportionality  of  activity  allows  for  algebraic 
operations  to  be  performed  on  this  set.  These  operations  allow 
for  the  interpretation  of  an  N dimensional  detector  response. 


Properties  of  Response  Sets 


Before  interpreting  a given  detector  response,  the  set  of 
responses  from  a specific  amount  of  activity  is  required  to  be 
constructed.  If  the  spatial  configuration  of  the  activity  is 
assumed  to  have  a specific  distribution,  such  as  a homogeneous 
distribution,  then  ignoring  statistical  uncertainty,  the 
response  set  will  be  composed  of  a point.  If,  however,  the 
distribution  is  not  known,  then  the  detector  response  results 
in  a number  of  points  where  each  point  on  the  response  set 
represents  one  or  more  spatial  distributions  of  activity. 

Although  the  location  of  the  source  within  the  container 
is  unknown,  the  source  may  be  restricted  to  a particular 
distribution.  For  example,  a radioactive  sample  vial  or  source 
planchet,  while  possessing  an  unknown  location  within  a room 
or  waste  container,  will  still  remain  in  its  original  shape. 
Therefore,  depending  on  the  restrictions  of  the  source  itself, 
different  response  sets  are  possible.  The  two  response  sets 
considered  in  this  dissertation  are  the  point  source  response 
set  and  complete  response  set.  The  point  source  response  set 
results  from  all  responses  from  a point  of  unit  amount  of 
radioactive  material.  The  complete  response  set  results  from 
all  possible  distributions  of  a unit  amount  of  radioactive 
material.  The  point  source  response  set  can  be  analytically 
constructed  by  forming  a parametric  set  of  detector  response 


equations  that  are  dependent  on  the  spatial  coordinates  of  the 
point  source.  For  N detectors  there  are  N equations  of  the 
form  f'~fn(x,y,z)  . The  complete  response  set,  however,  is  more 
difficult  to  construct  using  functional  relationships  since 
such  a construction  would  require  a systematic  search  of  the 
different  spatial  positions  possible.  Fortunately  as 
summarized  in  the  following  section,  the  complete  response  set 
can  be  derived  from  the  point  source  response  set.  In 
addition,  the  relation  of  the  complete  and  point  source 
response  sets  allows  for  the  complete  response  set  to  be 
considered  "between"  the  point  source  response  set. 

construction  of  the  Complete  Response  Set 

The  point  source  response  set  F is  the  set  of  all 
detector  responses  that  are  possible  from  one  unit  of  activity 
of  a single  point  source  and  is  defined  as 

F=  z)  for  all  xex,  yey  and  zez)  (2-1) 

where  X, Y and  Z represent  the  set  of  all  possible  locations  of 
a single  point  source.  If  this  unit  point  source  in  divided 
into  J amounts  a„a„ — , a,  at  respective  locations  (xt, y;/ 
the  detector  response  to  this  combination  of  subunit  sources 
is  defined  as 

a]f{x],y],z1) , where  £ Sj=l . 


(2-2) 


Due  to  the  formulations  that  result  from  the  above 
equation,  a careful  examination  of  Equation  2-2  is  in  order. 
Equation  2-2,  which  is  the  detector  response  resulting  from  a 
group  of  subunit  points  sources  adding  to  one,  is  equivalent 
to  a convex  combination  of  the  points  of  J points 

,f(.xI.yi.z2)  ,~.£{Xj.yj.  zj  . (2-3) 

This  equivalence  is  the  result  of  the  assumption  that  the 
detector  response  is  proportional  to  the  amount  of  activity 
measured.  From  Appendix  A,  the  set  of  all  convex  combinations 
of  the  elements  of  the  point  source  response  set  is  equivalent 
to  the  convex  hull  of  the  point  source  response  set.  Although 
the  examples  throughout  the  dissertation  present  a number  of 
simplified  formulations  for  /,  the  complexity  of  f and  its 
dependence  on  variables  such  as  location,  attenuation 
coefficient  of  the  matrix  and  scattered  radiation  does  not 
affect  the  above  relationship  as  long  as  the  response  remains 
proportional  to  activity. 

The  boundaries  of  the  complete  response  set  are  composed 
of  the  point  source  response  set  and  linear  combinations  of 
the  point  response  set.  For  a two  dimensional  response  space, 
the  boundaries  of  the  complete  response  set  will  be  the  point 
response  set  and  line  segments  whose  endpoints  are  elements  of 
the  point  response  set.  Example  2-1  formulates  a case  where 
the  extreme  points  of  the  point  source  response  set  define  the 
line  segment  boundary  of  the  complete  response  set.  These 
endpoints  represent  the  unit  source 


closest 


respective  detectors.  In  general,  an  N dimensional  complete 
response  set  will  have  the  boundaries  composed  of  the  point 
response  set  and  sections  of  N-l  dimensional  hyperplanes 
defined  by  elements  from  the  point  response  set.  The  question 
as  to  when  is  a response  "between"  other  responses  has  been 
resolved. 

If  the  number  of  quantities  and  respective  locations  of 
the  subunit  point  sources,  J,  approach  infinity,  then  the  set 
of  detector  responses  is  represented  as 

Ct-lcr:c,(x,y,z)  -Jjfr(x,  y,  z)  fix.  y,  zi  dxdydz  for  all  reR) 

where  R is  the  set  of  all  functions  integrable  on  the  space 
(X, Y, Z;  for  one  unit  of  activity.  The  above  equation  is  both 
a infinite  extension  to  the  definition  of  a convex  combination 
and  the  point  kernel  representation  to  a unit  amount  of 
radioactivity . 

The  advantage  of  demonstrating  that  the  complete  response 
set  is  the  convex  hull  of  the  point  source  response  set  is 
that  the  complete  response  set  can  be  formulated  in  terms  of 
the  point  source  response  set.  The  relationship  between  the 
sets  will  prove  advantageous  when  searching  over  the  complete 
response  set.  The  construction  of  a single  point  from  the 
complete  response  set  requires  N points  from  the  point  source 
response  set.  Each  of  the  N points  from  the  point  response  set 
has  a coefficient.  Since  the  coefficients  sum  to  one,  only  N-l 
coefficients  are  independent.  Therefore  2N-1  independent 


•(afe'TiS5,  lAj’Tipi1) 


of  t,  and  /, 


represented  in  equation  2-7.  The  area 


equation  of  this  line  is 
contained  by  the  point  source  response  set  and  the  line 
segment  is  the  complete  response  set.  Figure  2-4  illustrates 
the  complete  response  set. 

The  above  example  introduces  the  concept  of  linear 
equations,  at  least  partially,  bounding  a convex  set.  The 
hyperplane  that  bounds  a convex  set  C which  is  the  convex  hull 
of  a set  F , is  called  a supporting  plane.  A convex  set  will, 
at  each  extreme  point,  have  a supporting  plane.  Since  the 
intersection  of  convex  sets  is  a convex  set,  the  intersection 
of  the  supporting  plane  H with  the  convex  set  C or  C fl  H,  at 
the  point  x is  itself  a convex  set.  If  the  point  x is  a member 

combination  of  points  from  F.  These  points  from  F define  a 
hyperplane  P that  either  does  or  does  not  contain  the  set 
C 0 H.  If  P does  contain  the  set  C fl  H,  then  the  supporting 
plane  does  intersect  the  set  F.  If  P does  not  contain  the  set 
C n if,  then  some  direction  along  P will  form  a parallel 
hyperplane  defined  by  a point  from  F.  Therefore,  the 
hyperplane  H cuts  C and  does  not  bound  C and  H and  is  not  a 
supporting  plane.  The  supporting  plane  of  C will  contain  a 
point  from  the  set  F and  supporting  planes  of  a convex  hull  of 
a set  F will  intersect  the  set  F.  The  importance  of  the  above 
result  is  that  a linear  optimization  over  the  complete 
response  set  can  be  performed  over  the  point  source  response 
set  reducing  the  number  of  variable  required  for  a search. 


different  amounts  of  activity  and  different  source  locations. 


Figure  2-2.  Two  detector  system  distinguishing  between 
different  amounts  of  activity  and  different  source  locations. 


Figure  2-3.  Point  source  response  set  for  a line  source 
container. 


set  for  a line  source  container. 


Complete  response 


CHAPTER  3 

THE  LINE  INTERSECTION  (OR  RAV)  METHOD 
Introduction 

Previous  published  research  on  the  assay  of  randomly 
distributed  radioactive  material  attempted  to  fulfill  a 
lacking  in  the  subject  which  has  been  summarized  as  follows: 
“the  theoretical  ramifications  of  spatial  uncertainty  have 
never  been  analyzed  in  a systematic  way  nor  have  been  given  a 
unified  treatment.  Mathematical  techniques  for  rigorously 
assessing  the  spatial  uncertainty  of  a proposed  assay  system 
design  have,  until  recently,  been  lacking"  (Ben-Haim  1985b, 
pg.2). 

The  RMR  was  developed  as  a tool  for  assessing  the  spatial 
uncertainty  of  an  detector  assay  system.  As  acknowledged, 
however,  "careful  design  is  not  enough.  The  analyst  must 
properly  interpret  the  results  of  his  measurements"  (Ben-Haim 
1985b,  pg.5).  The  problems  of  assay  system  design  and 
interpretation  of  measurements,  however,  were  given  different 
mathematical  treatments  (Ben  Haim  1985b,  pg.  12) . While  the 
tools  of  assay  system  design  were  given  a deterministic 
treatment,  the  interpretation  of  detector  responses  was 
performed  utilizing  probabilistic  methods. 


The  reason  for  the  separation  is,  in  part,  historical. 
Articles  in  the  early  1980s  (Ben-Haim  and  Elisas  19B2,  Ben- 
Haim  1983,  Ben  Haim  et  al.  1981)  explored  the  use  of 
probabilistic  methods  to  interpret  detector  readings.  The 
probabilistic  interpretation  method  was  incorporated  with  the 
idea  of  constructing  sets  and  developing  the  concept  of  the 
RMR  (Ben-Haim  1983).  A number  of  articles  followed  that 
incorporated  the  idea  of  the  RMR  in  the  existing  set  theory 
(Ben-Haim  and  shanhav  1984,  Ben-Haim  1985a) . 

This  chapter  resolves  the  separation  of  design  and 
interpretation  by  proposing  a multi-detector  interpretation 
method.  Multi-dimensional  rays  in  conjunction  with  the 
response  set  are  utilized  for  the  purpose  of  developing  a 
method  that,  given  a set  of  detector  responses,  assesses  the 
amount  of  activity. 

Interring  a petector  ResB.onsg..from  & Known  Activity 

Throughout  the  dissertation  it  is  assumed  that  the 
detector  responses,  whether  from  a point  source  or  some  other 
distribution,  are  proportional  to  the  amount  of  radioactivity 
present.  Therefore,  if  a point  {£„£„...  ,£„)  is  the  response 
from  one  unit  of  activity,  then  {m£„m£„  . . . ,m£„)  is  the 
response  from  m units  of  activity. 

One  result  of  this  property  of  proportionality  is  that 
all  points  mf=(mf„m/3, . . .,mfw)  that  are  the  result  of  some  m 
units  of  activity  from  a specified  distribution  are  located  on 


the  same  ray  from  the  origin  with  a distance  from  the  origin 
proportional  to  the  m units  of  activity  present.  A point 
(•f;/  *1/  • • • / f«)  which  is  the  result  of  a unit  amount  of  source 
activity  is  located  on  the  ray  intersecting  the  origin  with 


d-Jf}*fj*. . . *f‘.  (3-2) 


system  response  is  (m The  slopes  of  this  ray 
are  identical  to  the  slopes  of  the  ray  from  one  unit  of 
activity.  As  shown  in  Figure  3-1  the  distance  from  the  origin 


d=J (inf,) 2* (mf2) J*. . . ♦ (mf„) 1 (3-3) 

or  a times  the  amount  from  one  unit  of  activity. 

interpreting  Activity  from  a Kntnm-Esg.Bsng.e_Po.lnt 


inferred.  The  problem  of  interpretation  is  to  calculate  the 
activity  that  results  in  a response  point  of  unknown  activity, 


but  identical  source  configuration,  utilizing  a response  point 
of  known  activity.  Since  the  points  of  known  activity  and 
unknown  activity  are  located  on  the  same  ray  and  the  distances 
from  the  origin  are  proportional  to  the  amount  of  activity 
present,  the  unknown  amount  of  activity  is  equal  to  the 
distance  from  the  point  response  of  the  unknown  activity 
divided  by  the  distance  from  the  response  point  of  the  known 
activity.  Given  a response  point  (g,,g., . . . ,g„)  from  m,  units 
of  activity  and  a response  point  from  an  unknown  amount  of 
activity  m,=  (h„Ji„  . . .,h„) , the  amount  of  activity  is 

.-<<**»■  ■■*«8  „„> 

A two-dimensional  example  is  shown  in  Pigure  3-2. 

Interpretation  of  Responses  from  a Unit  Response  Point 


For  a given  amount  of  activity,  a response  set  may 
possess  a range  of  possible  responses  along  a ray.  If  the 
activity  of  a given  detector  response  is  inferred  from  a range 
of  possible  responses,  the  interpretation  of  the  activity  from 
the  responses  is  performed  by  calculating  two  ratios.  These 
ratios,  in  both  cases,  have  as  the  numerator  the  distance  from 
the  origin  of  the  response  of  unknown  activity  (/)„&,,  . . .,JiK) . 
The  denominators  are  the  maximum  and  minimum  distances  from 
the  origin  to  the  responses  of  known  activity,  where  the 
are  taken  along  the  same  ray  th 


distances 


»l“^.  (3_6) 


C,-3lCl. 


where  for  g and  f,  VsHf,,  , ....  • 

The  point  source  response  set  can  be  formulated  as  the 
parametric  representation  of  the  spatial  coordinates  of  a unit 
point  source  with  a specified  source  container.  The 
interpretation  is  performed  over  the  complete  response  set 
which  is  the  set  of  all  convex  combinations  of  the  point 
source  response  set.  The  properties  of  the  point  source 
response  set  often  allow  the  complete  response  set  to  be 
constructed  with  less  complexity  than  resorting  to  calculating 
all  the  convex  combinations  in  the  point  source  set.  The 
following  examples  illustrate  the  application  of  the  line 
intersection  method  for  a number  of  simplified  detector  and 
source  configurations. 

Example  3-1.  The  radioactivity  is  confined  to  the  line 
segment  -O.SsxsO.5.  The  detectors  are  located  on  the  x axis 
at  x=l  and  x=-l  with  response  functions 

f,-  and  fa-  (i  1 )2  ■ (3-9) 

The  point  source  response  set,  which  is  determined  from  the 
above  two  response  functions,  is  illustrated  in  Figure  3-4. 
The  construction  of  the  complete  response  set  sometimes  can  be 
determined  from  the  behavior  of  the  point  response  functions. 
The  second  derivative  of  f2  with  respect  to  f,  is 


3--83S*- 


(f;/fj)  = (4,2) , the  slope  of  the  1 


-Has-- 


and  f?-0 . 5/, , or  (fr, t2) “(2.96,1.48) . As  demonstrated  in  Figure 
3-6,  the  maximum  and  minimum  activity  that  could  cause  the 
response  (4,2)  is  2.74  and  1.35,  respectively.  The  ratio  of 
the  maximum  activity  to  minimum  activity  is  2.027. 

The  ratio  of  the  maximum  to  minimum  activity  possible  is 
dependent  on  the  detector  response  that  is  initially  obtained. 
For  the  above  example,  if  the  given  detector  response  was 
(fJ,f,)  = (0.444,4) , the  ratio  of  maximum  to  minimum  activity  is 
one.  The  maximum  possible  value  of  the  ratio  is  2.22  which  is 
the  result  of  the  given  response  of  detector  #1  equal  to  the 
response  of  detector  #2.  This  maximized  ratio  is  equivalent  to 
the  RMR  discussed  previously.  The  utilization  of  the  line 
intersection  method  for  the  RMR  is  discussed  in  the  next 

The  above  example  can  be  extended  into  areas  such  as 
radioactive  waste  assessment.  An  applicable  problem  is  the 
assessment  of  activity  of  a 55  gallon  waste  container.  Waste 
container  assays  are  usually  performed  by  assuming  that  the 
activity  is  homogeneously  distributed  or  that  any 
heterogeneous  distributions  will  not  significantly  affect  the 
assessment.  If  it  is  no  longer  assumed  that  the  activity 
possesses  spatial  homogeneity,  the  ray  method  allows  for  the 
use  of  multiple  detector  measurements  to  reduce  the  assessment 
uncertainty. 

Example  3-2.  The  radioactivity  is  confined  to  an 
infinitesimally  thin 


disk  of  radius  R.  As  illustrated  in 


Figure  3-7,  two  detector  readings  are  performed  above  the 
centerline  of  this  disk.  The  problem  is  to  determine  the  range 
of  possible  activities  that  could  result  in  the  two  given 
detector  readings. 

Using  cylindrical  coordinates,  the  equation  of  the  disk 
that  confines  the  radiation  is  r=R.  The  detector  readings  are 
performed  at  z=Zj  and  z=Z3  with  . The  response  functions 


The  second  derivative  of  t , with  respect  to  /,  is 


d‘f3  _ -2 
dfi 

Therefore,  the  point  response  function  of  t , in  terms  of  t,  is 
concave  downwards  and  the  upper  boundary  of  the  complete 
response  set  is  the  point  response  set  and  the  lower  boundary 
of  the  complete  response  set  is  the  line  segment  connecting 
the  endpoints  of  the  point  response  set.  As  shown  in  Figure  3- 
8,  the  boundaries  of  the  complete  response  set  are  the  point 
response  set  and  the  line  segment 


I =(aft  (0)  ♦ (1-a)  f,  (1) , of,  (0)  * (l-o)  fj  (1) ) 
n • 11 -a)  ' _g_ , (I-*)  \ where  Osgsl. 


The  equation 


line  containing  this  segment  ii 


fj-wfj+h 

. 2,3(Z?*1) 


Therefore  f,  and  /,  in  terns  of  Z,  and  Z,  are 


f _ f _ h^Zl-2%)  (3-17) 

1 h;Zj  (Z2J *1)  -hjZ?(Z,J+l)  2 hjZ?  (Z*+l)  -AjZi  (Zi*l) 

A ray  from  the  origin  to  the  point  (h(/h3)  intersects  the  point 
source  response  set  at 


>h->h  f . A, -A, 

lh(zS-zl)  ' 3 h^z}-zl)  ' 


(3-18) 


As  a numerical  example,  the  two  detectors  readings  are  taken 
at  Z,=l  and  Zj=1.414.  The  detector  responses  are  given  as 
(h„hj)  = (l,0.55) . The  intersection  with  the  point  response  set 
is  (0.818,0.45).  As  demonstrated  in  Figure  3-9,  the 
intersection  with  the  line  segment  boundary  of  the  complete 
response  is  (f/,f2)“(0. 771, 0.424)  . The  range  of  activities  that 
could  result  from  the  given  detector  response  point  (1,0.55) 
is  1.222  to  1.297  units  of  activity,  a difference  of  6.1%.  If 
only  one  detector  was  available  then  the  range  of  activities 


possible  from  a detector  reading  of  hf"l  is  one  to  two  units 
of  activity,  a difference  of  100%. 

The  above  example  serves  as  a simplified  version  of  a 
number  of  applied  problems.  An  unknown  distribution  of 
radioactive  material  on  a sample  planchet  could  be  assessed  by 
counting  the  planchet  on  two  different  shelf  locations.  An 
aerial  survey  could  assess  the  amount  of  ground  contamination 
by  taking  measurements  at  two  different  altitudes.  Floor 
contamination  of  a room  may  be  assessed  by  survey  data  taken 
at  different  heights  above  the  floor. 

Example  3-3.  The  above  two  detector  problems  have  the 
luxury  of  the  use  of  graphical  technigues  for  visualization 
and  calculation.  For  greater  than  two  detectors  the  problems 
become  more  complex  due  to  increasing  number  of  constraints 
that  are  involved.  A three  detector  problem  utilizing  the  ray 
method  is  presented  in  the  following  example.  The  source  is 
confined  to  a cube  of  dimension  one  by  one  by  one  centered  on 
the  origin.  The  radioactivity  is  restricted  to  a single  point 
source.  Radiation  detectors  tl,  12  and  #3  are  located 
respectively  at  (1,0,0),  (-1,0,0)  and  (0,1,0).  The  three  point 
source  response  functions  in  terms  of  the  spatial  x,  y and  z 
coordinates  are  defined  as 

f i * — — r — - , f , ^ — - — - and  i , ■ — i — - — - . 

(l-x)**y1+z‘  (l*x)  ,*y**z2  xJ*(l -y)2*z2 

A given  line  passing  through  the  origin  is  f,=s,fj  and 


£■■  VW..  A W,.(3-24) 


^ ffj  (xvy.  Z)  s m* smax»  ?I  (x?y,  z)  ■ 


7lU,y,X)  ' <!-*> ’**’*’-*>  (Li,’)  *2 


-<z*+l>  -<**♦!)  j.(’  28> 


As  a second  numerical  example,  the  detector  responses  are 
gives  as  (h„h2,h2)  = (2. 909, 1.185,2. 909)  . The  constants  for  the 
problem  are  s,=2.4SS,  Sj=l,  b-l  and  w=2.375.  The  applicable 
range  of  x is  based  on  the  following  equation 

X- 

which  allows  x to  possess  a range  of  .233  to  .301.  The  amount 
of  activity  present  is 

l.SSSsia^.dOS.  (3-30) 

In  summary,  the  proposed  N detector  interpretation 
problem  has  an  N vector  response  point  and  an  N dimensional 
unit  response  set.  The  response  point  defines  a ray  from  the 
origin  composed  of  N-l  linear  equations,  which  intersects  the 
N dimensional  response  set.  The  distance  of  the  given  response 
point  from  the  origin  divided  by  the  distance  from  the  origin 
of  the  endpoints  of  the  line  segment  intersecting  the  unit 
response  set  give  the  maximum  and  minimum  amounts  of  activity 
that  could  result  in  the  given  response  set.  Due  to  the  linear 
relationship  between  the  detector  responses,  the  response  of 
only  one  detector  is  required  for  the  calculation  of  the 
ratios.  The  N-l  linear  equations  allow  for  the  formation  of  N- 
1 relationships  between  the  spatial  variables  (x,y,z).  These 
linear  equations  limit  the  permissible  ranges  of  the  spatial 
variables.  In  addition,  these  linear  equations  can  reduce  the 


number  of  dependent  spatial  variables  of  the  detector  response 
function. 

In  order  to  determine  the  range  of  possible  activity,  a 
search  procedure  must  determine  the  maximum  and  minimum 
detector  response,  subject  to  the  constraints  of  the  problem. 
In  terms  of  the  spatial  variables,  both  the  detector  responses 
and  the  constraints  are  nonlinear  functions.  Ho  general  method 
exists  for  optimizing  a nonlinear  function  subject  to 
nonlinear  constraints.  Utilizing  a local  search  procedure  for 
the  above  type  problem  may  result  in  a local  optimum  that  is 
not  the  global  optimum  of  the  problem.  In  terms  of  detector 
space  however,  the  search  of  the  interpretation  problem  is 
constrained  to  a line  segment.  The  endpoints  of  the  line 
segment  contain  the  maximized  and  minimized  detector 
responses . 

The  interpretation  problem  is  defined  ultimately  in  terms 
of  the  spatial  variables  which  do  not  linearly  constrain  the 
problem.  In  order  to  utilize  the  linear  constraint  of  the 
problem,  the  following  search  procedure  is  recommended.  After 
the  detector  response  appears  to  be  optimized  in  terms  of  the 
spatial  variables,  the  optimized  detector  response  is 
incremented  in  the  direction  (with  respect  to  detector  space) 
that  the  search  has  previously  progressed.  The  spatial 
variables  that  are  calculated  from  the  incremented  detector 
response  are  determined  to  be  within  or  without  the 
permissible  ranges.  If  the  variables  satisfy  the  ranges,  then 


the  calculated  response  is  not  an  optimum  solution.  The 
incremented  response  is  used  as  an  initial  point  £or 
restarting  the  search. 

In  order  to  include  unconstrained  distributions,  the 
search  is  required  to  be  performed  over  the  complete  response 
set.  The  search  over  the  complete  response  set  introduces 
additional  variables.  Examples  3-2  and  3-3  constructed  the 
boundaries  of  the  complete  response  set  by  determining  the 
curvature  of  the  point  source  response  function.  If  the 
curvature  were  solely  positive  or  negative,  an  upper  of  lower 
boundary  is  determined.  A more  general  method  of  searching 
over  the  complete  response  set  is  performed  by  generating  a 
detector  response  point  from  the  complete  response  set  from  N 
points  of  the  point  response  set  in  conjunction  with  N- 1 
independent  coefficients. 

A set  of  detector  responses  from  a two  detector  point 
source  response  set  is  represented  as  (/,,/j).  A point  from  the 
complete  response  set  is  represented  as  (c,,c2)  or  in  terms  of 
(/„/,)  as 

c, fl  and  c,*Bfj+(l-#)fa.  (3-31) 

If  a given  detector  response  point  (h„h3)  results  in  the  slope 
s then 

c2-ofs+(l-«)  £j-OSfj+(l-B>  sf,'.  (3-32) 

In  terms  of  s and  the  points  (f„f ,)  and  from  the 

point  source  response  set,  alpha  is 


sfl-fl 


uSj*-fSS-  ”' 


a = L-Li-y.) 1 d-*')2 (3-38) 

(!-*')  2 ' <1**)>  + (l**1)2  " d-x') 

The  search  over  x and  x'  is  reduced  to  a search  over  x by 
recognizing  that  if  the  complete  response  set  can  be  generated 
from  a fixed  extreme  point  on  the  point  response  set,  such  as 
(/,,f2)  = (. 444,4),  then  a variable  point  on  the  point  source 
response  set  and  the  line  segment  connecting  the  two  points  is 
sufficient  to  define  the  complete  response  set.  If  the 
variable  point  is  (f,,f3)*(4,  .444) , the  line  segment  connecting 
this  point  with  the  fixed  point  (.444,4)  creates  part  of  the 
boundary  of  the  complete  response  set. 

If  the  x'-.5,  then  alpha  is 


where  Osasi . 


An  alpha  outside  the  range  of  (0,1)  indicates  that  the 
intersection  of  the  line  created  from  the  variable  and  fixed 
response  points  and  the  ray  created  from  the  given  response 
point  of  unknown  activity  occurs  outside  the  line  segment  and 
therefore  outside  the  complete  response  set. 

If  the  detector  response  from  an  unknown  amount  of 
activity  is  given  as  (h(/hj)-(2,4)  (s=2) , alpha  is 


The  range  of  the  amount  of  activity  that  could  result  in  the 
detector  response  (2,4)  is  written  as  the  following 
constrained  optimization 


(!-*>»'  (!♦*)* 


Figure  3-10  illustrates  the  search  for  a fixed  point  (0.444,4) 
and  a variable  point  (1,1)  resulting  in  the  point  from  the 
complete  response  set  (0.67,1.72). 

The  constraint  limits  the  range  of  x to  between  x=-0.5 
(alpha*. 29)  and  x=. 17  (alpha=l) . The  amount  of  activity  that 
could  be  present  is 

■ 29,  71  i I (3-42) 

.5*  1.5’J 


1.35SU!4S2.74  (3-43) 
giving  the  same  result  as  found  in  example  3-1. 

If  the  unit  set  of  detector  responses  are  generated  from 
a computer  code,  the  search  procedure  can  not  take  advantage 
of  any  analytical  relationships.  Without  the  advantage  of 
analytical  relationships  to  simplify  the  search,  the  assessor 


is  confronted  with  the  choice  of  generating  elements  of  the 
complete  response  set  from  varying  spatial  distributions  or 
generating  only  elements  of  the  point  source  response  set.  If 
only  the  point  source  response  set  is  generated,  then  the 
search  over  the  complete  response  set  is  performed  by  either 
convex  combinations  of  the  point  source  response  set  or 
evaluating  the  linear  N-l  dimensional  hyperplanes  which  define 
the  boundaries  of  the  complete  response  set.  If  the  boundaries 
of  the  complete  response  set  are  determined,  the  search  is 
simplified  to  the  point  response  set  and  the  linear 
boundaries.  If  the  search  is  performed  by  convex  combinations, 
then  the  number  of  independent  variables  in  the  search  is 
equal  to  the  number  of  detectors  present.  Although  the  size  of 
the  generated  set  determines  partially  the  accuracy  and  CPU 
time  necessary  to  optimize  the  responses,  the  same 
consideration  illustrated  in  the  evaluation  of  the  above 
simplified  examples,  that  of  exploiting  the  properties  of  the 
point  response  set  in  order  to  simplify  the  search  of  the 
complete  response  set,  is  an  excellent  tool  for  the  reduction 
of  the  calculational  effort. 
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Figure  3-1.  Relationship  between  a response  from  a unit 
activity  and  a response  from  a unit  activity  multiplied  by  m. 
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Figure  3-2.  The  amount  of  activity  from  point  h that  is 
inferred  from  point  g is  ms(dh/dgj . 
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Figure  3-3.  The  amount  of  activity  from  point  k that  is 
inferred  from  the  unit  activity  line  segment  which  connects 
points  g and  h is  Sm,S  . 


Figure  3-4.  Point  source  response  set  for  a line  source 
container. 


a line  source  container. 


Complete  response  set  for 


Figure  3-6.  Line  intersection  method  with  a two  detector 
complete  response  set. 


configuratic 


Figure  3-8.  Complete  response  set  for  two  detectors  with  disk 


Figure  3-9.  The  line  intersection  method  for  the  complete 
response  set  of  two  detectors  with  a disk  source  container. 


Figure  3-10.  Fixed  point  search  procedure  with  variable  point 
at  (1,1)  and  fixed  point  at  (4,. 444). 


CHAPTER  4 

RELATIVE  MASS  RESOLUTION 

introduction 

The  previous  chapter  expanded  the  application  of 
dimensional  space  to  include  the  interpretation  of  multiple 

detector  responses  result  in  a different  ranges  of  interpreted 
activity.  To  quantify  this  range,  the  ratio  of  the  maximum  to 
minimum  amount  of  activity  of  this  range  is  calculated. 

The  application  of  dimensional  space  has  been  previously 
applied  in  published  work  to  the  problem  assessing  the  maximum 
uncertainty  of  a detector  system.  This  uncertainty  is 
determined  by  the  calculation  of  the  RMR  which  has  been 
defined  as  the  ratio  Cm* A)/m  where  (m*A)  is  the  smallest 
quantity  of  radioactive  material  greater  than  an  amount  m of 


a performance  criterion  could  assist  both  the  analyst  and  the 
designer  detector  systems. 

A previously  developed  algorithm  (Ben-Haim  1983),  which 
utilizes  supporting  planes,  determines  the  RMR  by  searching 
over  the  point  source  response  set  versus  the  complete 
response  set,  thus  reducing  the  number  of  variables  required 
for  the  search.  The  method  is  restricted  to  the  case  where  the 
response  set  is  convex.  Since  the  complete  response  set  is  the 
convex  hull  of  the  point  source  response  set,  the  restriction 
does  not  pose  a barrier  to  the  solution  of  the  RMR  for 
unconstrained  source  distributions.  If  the  possible  set  of 
source  distributions  is  restricted,  however,  then  the 
resulting  response  set  might  not  be  convex.  A source 
restricted  to  a single  point  source  results  in  a response  set 
which,  in  general,  is  not  a convex  set.  In  addition,  the  above 
search  procedure  does  not  guarantee  that  a local  optimum 
solution  is  a global  optimum  solution. 

In  order  to  address  the  above  shortcomings,  the  line 
intersection  method  of  the  previous  chapter  is  extended  to 
include  a method  for  searching  for  the  RMR.  The  ray  method,  in 
addition  to  incorporating  the  same  formulations  that  were  used 
in  interpreting  a detector  response,  is  able  to  search  for  the 
RMR  over  nonconvex  response  sets  and  guarantee  that  a local 
optimum  is  a global  optimum  for  convex  response  sets. 

The  number  of  independent  variables  required  for  the 
search  of  the  RMR  is  dependent  on  the  number  of  measurements. 


If  relationships  can  be  developed  between  the  slopes  of  the 
upper  and  lower  supporting  planes  and  the  response  set,  the 
number  of  independent  variables  can  be  reduced.  A method  of 
utilizing  the  derivatives  of  the  point  source  response  set  has 
been  developed  which  can  reduce  both  the  number  of  variable 
and  the  number  of  required  optimization  steps. 

Unlike  the  problems  presented  in  this  dissertation,  the 
calculation  of  most  multiple  detector  problems  will  require 
computer  based  search  methods.  If  local  type  search  methods 
are  utilized,  local  optimum  are  not  necessarily  the  global 
optima  of  the  solution.  Therefore  the  development  of  methods 
that  guarantee  the  global  applicability  of  local  solutions 
have  been  investigated.  The  line  intersection  method  for  the 
RMR  is  demonstrated  to  be  a global  optimization  method  for 
convex  sets.  Finally,  a technique  has  been  introduced  into  the 
supporting  plane  search  procedure  to  allow  for  the  global 
optimization  of  convex  sets. 

Basic  Formulation  of  the  RMR 
Distinguishing  Activity  Along  a Ray 

The  description  of  the  RMR  has  been  previously  formulated 
in  terms  of  N-l  dimensional  hyperplanes.  The  RMR,  however,  can 
be  described  in  terms  of  distances  from  the  origin.  The 
following  discussion  of  the  RMR  dovetails  two  fundamental 
descriptions  of  the  RMR:  the  N-l  dimensional  hyperplane 
approach  from  previous  literature  and  the  distance  from  the 


origin  approach  developed  in  the  previous  chapter.  The 
addition  of  the  distance  to  the  origin  approach  allows  for  the 
solution  of  problems  which  the  method  of  hyperplanes  does  not 
address.  Before  discussing  multi-detector  formulations  and  the 
techniques  developed,  the  following  single  detector  example 
problem  illustrates  that,  although  the  concept  of  the  RMR 
appears  foreign  to  practitioners  of  waste  assay,  the  RMR,  at 
least  for  the  one  dimensional  response  case,  is  an  intuitive 
concept  for  measuring  the  performance  of  a detector 
configuration. 

Example  4-1.  The  radioactivity  is  confined  to  the  line 
segment  -.Ssxs.S.  The  detector  is  located  on  the  x axis  at  x=l 
with  response  function 


The  detector  response  to  a point  source  of  unit  activity  is 
from  0.44  to  4 inclusive.  If  the  unit  point  source  is  divided 
into  a number  of  subunit  point  sources,  any  possible  detector 
response  resulting  from  a distribution  of  these  subunit 
sources  will  be  between  0.44  and  4.  Assuming  proportionality 
of  detector  response  to  activity,  nine  units  of  activity  will 
result  in  a response  between  4 and  36.  If  the  detector 
response  for  some  unknown  amount  and  distribution  of  activity 
is  given  as  four,  then  the  possible  amount  of  activity  present 
is  between  one  unit  and  nine  units.  Therefore,  the  RMR  of  this 
single  detector  system  is  nine  over  one  or  nine.  However, 


instead  of  calculating  the  RMR  from  its  formal  definition 
( the  RMR  can  be  determined  by  calculating  the  ratio 
of  the  maximum  and  minimum  responses  from  a set  of  responses 
with  a fixed  activity.  In  the  above  example,  the  RMR  can  be 
calculated  from  only  the  unit  response  set.  The  maximum  over 
the  minimum  response  for  one  unit  is  4/0.44  or  9. 

Since  a single  unit  source  will  result  in  more  extreme 
single  detector  responses  than  a divided  unit  source,  the 
ratio  of  the  maximum  to  minimum  detector  response  from  a point 
source,  at  least  for  a single  detector  response  with  the 
assumption  of  proportionality  of  response  to  activity, 
determines  the  RMR.  Part  of  the  problem  in  extending  the  above 
example  to  multiple  measurements  is  the  confusion  as  to  the 
representation  of  an  extreme  response  and  how  point  source 
distributions  can  "contain"  other  distributions. 

The  application  of  dimensional  space  for  multiple 
detector  responses  allows  the  definition  of  a maximum  and  a 
minimum  detector  response  to  be  extended  to  multiple  detector 
configurations.  The  criterion  that  determines  the  magnitude  of 
a multiple  detector  response  is  the  distance  of  a multiple 
detector  response  from  the  origin  in  response  space.  The 
application  of  the  complete  response  and  point  source  response 
set  allows  for  the  relationship  between  distributions  of 
subunit  sources  and  unit  point  sources  to  be  extended  to 
multiple  detector  configurations. 


The  line  intersection  method  of  the  previous  chapter 
considered  a line  segment  resulting  from  a ray  intersecting 
the  complete  response  set.  Such  a segment  represents  different 
detector  responses  due  to  various  spatial  distributions  of  a 
fixed  amount  of  activity.  In  addition,  the  line  segment 
represents  proportional  detector  responses  which  are  the 
result  of  fixed  spatial  distributions  of  various  proportional 
amounts  of  activity.  For  example,  a detector  response  twice 
the  distance  from  the  origin  of  another  detector  response  is 
the  result  of  twice  as  much  activity.  Therefore,  the  set  of 
detector  responses  along  the  line  segment  are  the  result  of 
two  different  means:  different  distributions  with  fixed 
amounts  of  activity  and  fixed  distributions  with  different 
amounts  of  activity. 

The  calculation  of  the  RMR  along  a segment  of  a ray  from 
the  origin  can  be  considered  a one  dimensional  detector 
problem.  The  distances  of  the  endpoints  from  the  origin  of  the 
line  segment  are  the  maximum  and  minimum  responses.  The  RMR  is 
calculated  by  dividing  the  distance  of  the  endpoint  that  is 
further  from  the  origin  by  the  distance  that  is  closer  to  the 
origin.  In  terms  of  distances  from  the  origin,  the  equation  of 


As  mentioned  before,  supporting  planes  can  be  utilized  to 
determine  the  RMR.  A result  of  the  property  of  proportionality 
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the  resulting  equation  for  the  RMR  of  either  the  supporting 
plane  method  or  the  distance  from  the  origin  method, 
simplifies  to 


Distinguishing  Activity  within  a Set 

The  calculation  of  the  RMR  is  usually  not  confined  to  a 
single  segment  of  a ray  from  the  origin  but  to  a set  of 
responses,  such  as  the  complete  response  set.  The  complexity 
of  the  search  for  the  RMR  over  a set  versus  over  a ray 
increases,  since  the  RMR  calculated  for  one  particular  line 
segment  may  be  different  than  the  RMR  calculated  for  another 
line  segment.  If  a response  set  is  multiplied  by  a ratio  that 

contain  points  in  common  with  the  initial  set.  The  RMR  can  be 
calculated  by  searching  for  the  line  segment  which  results  in 
the  greatest  maximum  to  minimum  ratio.  If  the  calculation  of 
the  RMR  is  formulated  in  terms  of  distances  from  the  origin, 
the  procedure  of  finding  the  RMR  over  a complete  set  is 
similar  to  the  ray  method  outlined  in  Chapter  3 with  the 
addition  that  the  calculation  is  performed  over  all  possible 
rays  from  the  origin  that  intersect  the  response  set,  C.  The 
ray  method  formulation  of  the  RMR,  over  a set  C,  is 
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If  the  response  set  is  convex,  the  method  of  supporting 
hyperplanes  can  be  utilized  to  determine  the  RMR.  The 
distances  from  the  origin  of  the  upper  and  lower  supporting 
planes  for  given  vector  a=(a„aj. . .,aK)  and  a set  C are 
calculated,  as  illustrated  in  figure  4-2.  The  ratio  of  the 
distances  of  these  parallel  planes  serves  as  the  trial  factor, 
mT,  of  the  RMR.  The  set  calculated  when  each  point  in  C is 
multiplied  by  the  trial  factor  mT  is  mfi.  The  response  points 
ntjC  are  located  on  upper  supporting  plane  of  C and  the  lower 
supporting  plane  of  mfi.  The  sets  C and  m,c  are  separated  by 
the  upper  supporting  plane  of  C (or  the  lower  supporting  plane 
of  mfi) . The  only  location  for  identical  response  points 
between  the  sets  C and  m-fi  is  on  the  upper  supporting  plane  of 
the  set  C.  If  a point  from  the  set  C is  multiplied  by  mr,  the 
point  nijC  is  located  on  or  above  the  upper  supporting  plane. 
If  the  sets  mfi  and  C do  not  have  points  in  common  along  the 
upper  supporting  plane,  then  the  RMR  is  smaller  than  the  trial 
factor.  Due  to  the  property  that  convex  sets  contain  all 
convex  combinations  of  their  elements,  the  smallest  trial 
factor  from  the  different  possible  slopes  of  the  supporting 
planes  is  the  RMR  (Ben-Haim  1985a) . The  equation  for  the  RMR 
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The  correct  slope  of  the  supporting  planes  of  e 
set  C is  illustrated  in  Figure  4-3.  If  the 
convex,  then  the  above  equation  does  guarantee  t! 


If  the  response  set  C is  the  convex  hull  of  a set  F,  the 
search  for  the  supporting  planes  can  be  performed  over  the  set 
F versus  the  set  C.  The  elements  g and  h in  the  above  equation 
are  then  required  to  belong  to  F instead  of  c.  since  the  point 
source  response  set  is  easier  to  formulate  than  the  complete 
response  set,  the  method  of  supporting  planes  is  less  complex 
than  the  method  utilizing  rays  intersecting  the  response  set. 

An  algorithm  based  on  the  utilization  of  supporting 
planes  for  the  evaluation  of  the  RMR  has  been  considered  as  an 
inner  and  outer  search  (Ben-Haim  1983) . The  inner  search 
locates  the  upper  and  lower  supporting  planes  for  a given 
(a„a„ — ,a„)  by  maximizing  and  minimizing  the  numerator  and 
denominator  of  the  above  equation.  The  upper  supporting  plane 
of  the  set  F is  also  the  upper  supporting  plane  of  the  set  C, 
while  mint<a,gr>  is  the  lower  supporting  plane  of  the  set  C. 
The  minimization  of  the  ratio,  maxa<a/h>/minx<a,g>,  over  a is 
the  "outer"  search.  The  number  of  independent  variable  in  the 
total  search  is  3N-1  where  N is  the  number  of  detectors 
present . 


A number  of  difficulties  are  present  in  the  application 
of  the  algorithm.  While  a gradient  type  optimization  has  been 
recommended  as  a search  method,  both  the  "inner"  and  "outer" 
searches  suffer  from  the  problem  that  a local  optimum  is  not 
necessarily  of  global  optimum  (Scott  1991) . In  addition,  in 
order  to  prevent  the  algorithm  from  resulting  in  solution  of 
a set  of  planes  that  "straddle"  the  origin,  a rejection 
criterion  is  established  that  an  optimum  solution  must  have 
agreement  between  the  signs  of  the  numerator  and  denominator 
(Ben-Haim  19S3,  pg.  58)  . No  follow-up  procedure  is  given, 
however,  if  an  optimized  solution  must  be  discarded  for 
nonagreement  of  signs  (Scott  1991) . In  addition,  the  search, 
if  implemented  as  described,  could  result  in  zero  denominators 
(Scott  1991) . 

The  supporting  plane  search  method  for  the  RMR,  although 
not  a global  optimization  technique,  possesses  the  advantage 
of  having  the  search  space  limited  to  the  point  source 
response  set  (Ben-Haim  1985b) . The  search  is  performed  over  a 
smaller  search  space  than  the  ray  method.  The  result  of  the 
smaller  search  space  is  that  the  maximization  and  minimization 
of  the  respective  numerator  and  denominator  is  performed  over 
only  one  point  and  W-l  weight  coefficients  for  the  supporting 
plane  method  versus  N points,  N-l  weight  coefficients  and  the 
slope  of  the  rays  that  intersect  the  response  set. 


Use  of  Derivatives  as  Weight  Coefficients 


If  a mathematical  relation  between  the  "inner"  and 
"outer"  search  is  constructed,  the  number  of  independent 
variables  and  the  number  of  optimization  steps  can  be  reduced. 
If  the  point  source  response  function  is  a convex  function, 
all  linear  combinations  of  the  point  source  response  set  will 
be  greater  than  or  equal  to  the  point  source  response  set.  If 
the  derivative  exists  where  the  lower  supporting  plane 
intersects  the  response  set,  the  slope  coefficients  of  a lower 
supporting  plane  can  be  written  in  terms  of  the  derivatives  of 
the  tangent  point.  Since  the  upper  supporting  plane  is 
parallel  to  the  lower  supporting  plane,  the  slopes  upper 
supporting  plane  can  be  defined  in  terms  of  the  derivatives  of 
the  point  source  response  set.  starting  from  the  equation  for 


RMR=min.(°>aX|l(ai^ll*a2^2*’  ’ ' (4-55) 
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if  the  numerator  and  denominator  of  the  above  equation  is 
divided  by  a„,  the  resulting  equation  is 


The  directional  derivatives,  a_/a„  of  the  plane  can  be  written 
in  terms  of  the  derivatives  of  the  tangent  point  as 
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The  value  of  r!  at  the  above  constraint  is 


Substituting  .366  for  z2,  results  in  {f„/j)=(l,  .5)  and  an  RMR 

Utilizing  derivative  information  for  cases  where  the 
point  source  response  set  is  convex  or  concave  allows  for  an 
elimination  of  the  slope  coefficients  of  both  the  supporting 
hyperplanes  and  the  "inner"  search.  If  the  point  source 
response  set  is  not  strictly  convex  or  concave  the 
mathematical  relationships  formed  from  an  investigation  of  the 
derivatives  could  still  allow  for  the  a reduction  in  the 
number  of  local  search  variables  required. 

Use  of  the  Rav  Method  for  Nonconvex  Sets 

The  supporting  plane  method  only  searches  for  the  RMR  for 
a convex  set  (Ben-Haim  1983) . If  the  set  of  spatial 
distributions  is  constrained,  the  resulting  response  set  may 
not  be  convex.  A common  constraint  is  to  restrict  the  activity 
to  a single  point  source.  The  formulation  for  the  RMR  over  a 
point  source  response  set  is 
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For  any  r between  zero  and  one,  a value  for  cos8  can  be  found 
that  satisfies  the  above  constraint.  In  addition,  for  r=l, 
cos8  is  indeterminate.  Therefore,  ignoring  the  constraint,  the 
above  equation  for  the  RMR  is  maximized  when  r,=l  and  r;=Q , 
resulting  in  an  RMR  equal  to  two.  When  r= 0,  (f„f,)»(l,  .25)  and 
the  slope  s of  the  ray  that  determines  the  RMR  is  0.25.  Figure 
4-5  illustrates  the  point  source  response  set  and  the  ray  used 
for  calculating  the  RMR. 

Example  4-6.  The  source  is  confined  to  a cube  centered  on 
the  origin  with  sides  equal  to  one.  The  radioactivity  is 
restricted  to  a single  point  source.  Radiation  detectors  #1,#2 
and  #3  are  located  respectively  at  (1,0,0),  (-1,0,0)  and 
(0,1,0).  The  three  point  response  functions  in  terms  of  the 
spatial  x,  y and  z coordinates  are 

1 (1  -x)a*y‘*z3'  (l*x)**y**z*  a"d  xJ»(l-y)  J*zJ  ' 

A given  line  passing  through  the  origin  is  f,=s,f * and 
fi=sA  or  in  terms  of  the  spatial  coordinates 
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where  w-  **** , b-  2S^*3**1 , -.Ssy,z,s.5.  (4-101) 
The  solution  to  the  above  equation,  which  will  not  be  pursued 
further,  would  require  more  computational  effort  than  the 
other  example  problems  that  have  been  presented. 

Techniques  for  Global  Optimization  of  the  RMR 

The  search  for  the  RMR  utilizing  the  above  method  of 
substitution  of  derivatives  for  the  weight  coefficients 
assists  in  forming  a relationship  between  the  inner  and  outer 
search;  however,  the  problem  that  the  gradient  search  method 
locates  a local  optimal  solution,  but  not  a global  optimum 
solution,  still  exists.  Assume  that  a hyperplane  is  locally 
optimized  on  the  point  source  response  set.  If  the  hyperplane 
does  not  intersect  the  point  source  response  set  at  any  other 
location  then  the  hyperplane  is  globally  optimized  and  is  a 
supporting  plane.  If  however  the  hyperplane  does  intersect  the 
point  source  response  set  then  the  hyperplane  has  not  been 
globally  optimized.  The  additional  point  or  points  of 
intersections  of  the  nonoptimized  hyperplane  will  assist  in 
determining  the  global  optimization  of  the  hyperplane.  If  the 
hyperplane  intersects  the  point  source  response  set  then 

a,f1*ajf2*. . . *a0f„=b  (4-102) 
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for  the  points  that  satisfy  the  hyperplane  equation  with  a 
small  given  error.  The  point  that  satisfies  the  hyperplane 
equation  with  the  least  error  would  form  the  new  starting 

The  disadvantage  of  the  ray  method  compared  to  the 
supporting  plane  method  for  calculating  the  RMR  of  a convex 
set  is  that  the  ray  method  must  search  over  the  complete 
response  set  while  the  supporting  plane  method  searches  over 
the  point  source  response  set.  In  addition,  the  ray  method  is 
a constrained  optimization  while  the  supporting  plane  method 
has  been  presented  as  an  unconstrained  problem.  As  the 
following  discussion  illustrates,  if  the  line  intersection 
method  is  applied  to  a convex  set,  a locally  optimized  maximum 
to  minimum  ratio  is  a global  optimum. 

Consider  a ray  originating  from  the  origin  intersecting 
a convex  set.  This  intersection  results  in  a line  segment  of 
length  L and  a point  that  is  a maximum  distance  D.lul  from 
the  origin  and  a point  P„IB  that  is  a minimum  distance  from 
the  origin.  The  ratio  of  the  maximum  and  minimum  distances 
is  the  baseline  for  the  RMR  search.  If  the  ratio  is  a 
local  maximum  then  the  line  segments  in  the  neighborhood  of 
the  ray  of  the  local  optimum  have  smaller  ratios  than  the 

initial  ratio  D / D— . . If  it  is  assumed  that  the  calculated 

ratio  is  not  the  global  maximum  ratio,  then  there  exists 
another  line  segment  with  points  and  P'mi>  and  a ratio 

D'^/D'n*,  that  is  greater  than  the  ratio  Dm/D„.  Therefore, 


the  line  segment  connecting  P'^,  and  P'm  is  greater  than  the 
line  segment  connecting  Pm  and  PM.  since  C is  a convex  set, 
the  area  defined  by  the  linear  combination  of  points  P_,  p*,, 
P'nu,  and  P ' M is  contained  within  the  set  c.  The  length  of  line 
segments  contained  in  the  convex  sets  from  rays  originating 
from  the  origin  increases  from  where  the  ray  intersects  P^,  to 
where  the  ray  intersects  P'm.  However,  it  was  assumed  that 
the  lines  segments  in  the  neighborhood  of  the  line  segment  P^ 
to  P^  are  smaller  than  the  above  line  segment.  Therefore,  the 
assumption  that  a global  optimum  exists  that  is  different  than 
the  local  optimum  is  incorrect. 

As  stated  previously,  the  ray  method  searches  over  the 
convex  hull  instead  of  only  the  point  source  response  set.  If 
. . .,1K)  are  members  of  the  point  source  response  set 
F,  the  convex  hull  is  defined  as 


c-£  «,f,) , where  £ a, 
onvex  set  in  terms  of 


l and  a al.  I 
he  ray  method 


where  for  c a 
After  substituting  ■ 


i constraints  1 


above  equation 


Ca-s^i  and  c„-^  aBf‘ 


— (*)• 


and  Cj-aVd-a)/,'.  (4-108) 

response  point  (h„h,)  results  in  the  slope 


Os/,'. 

nd  i 


(4-110) 


(4-111) 


c’’  *(  u^)' 


(l«x')»"  tl-x')a 


71^ '1T^  'TT^T 


the  two  points.  If  the  variable  point  is  (f„/j)  = (4, 0.444) , the 
line  segment  connecting  this  point  with  the  fixed  point 
( .444,4)  creates  part  of  the  border  of  the  complete  response 


f 3-1)  *2x(s*ll  -x3 


5S-2.25) 


9-1)  .. 


<1*X)*(1-X>*  '*  >‘-S-2.25)  (4-117) 

where  Osasl. 

An  alpha  outside  the  range  of  (0,1)  would  indicate  that  the 
intersection  of  the  line  defined  by  the  variable  and  fixed 
response  points  and  the  ray  created  from  the  given  response 
point  of  unknown  activity  occurs  outside  both  the  line  segment 
and  the  complete  response  set. 

If  the  initial  starting  point  of  s is  given  as  s=2,  alpha 


• where  osasl. 


(l*x)*(l-x)* 

The  inequality  constraint  limits  the  range  of  x to  between 
x*”.5  where  alpha=0.29,  (f,,/,)  = (4,0. 44)  and  x-0.17  where 

alpha=l  and  731,1.45).  The  equation  for  the  RMR  is 


RMR-max.max.  f whete  Ososl.  (4-119) 

Substituting  the  equations  for  alpha,  t,  and  g,  gives 


(l*x>*(l-x)*" 
where  Osasl. 


For  s=2  the  above  equation  is  maximized  for  x=-0 . 5 , [f,,f2)  = 
(4,. 44),  and  x'=0.17  (f„f2)=(0.73,1.45) . As  illustrated  in 
Figure  4-7,  the  optimized  detector  response  of  the  numerator 
is  the  extreme  point  of  both  the  point  source  response  set  and 
the  source  container.  The  optimized  response  of  the 
denominator  has  an  alpha  equal  to  one  indicating  that  the 
minimized  response  point  from  the  complete  response  set  is  on 
the  point  source  response  set.  The  ray  intersects  both  the 
line  segment  that  connects  the  two  extreme  points  of  the  point 
response  set  and  the  point  source  response  set  itself. 
Possibly,  the  optimized  ray  will  intersect  the  above  line 
segment  and  the  point  source  response  set.  Therefore,  the 
numerator  is  kept  constant  at  x=-0.5  while  the  denominator  is 
allowed  to  vary.  The  optimization  over  the  complete  response 
set  of  the  denominator  is  changed  to  an  optimization  over  the 
point  source  set.  The  resulting  equation  is 

where  Osssl. <4-121) 

The  intersection  of  the  point  source  response  set  with  a ray 
defined  by  determines  the  following  x. 


(s-l)x2-2(s+l)x*(s-l) 


-4-«9S- 


J 2«Wl”  -?.‘*2  (/»•)■  “-1”1 


(Cj.e,)  =(«f1(0)  + U-«)fl(l)/«f!(0)t(l-«)/s(l)) 


■-*&— ^ssr 


point  (h„Jij)  intersects  the  boundaries  of  the  set  C on  the 
point  source  response  set  and  the  line  segment  connecting  the 
extremes  of  the  point  source  response  set.  The  ratio  of  the 
distances  from  the  origin  of  the  intersection  of  the  line 
segment  with  the  intersection  of  the  point  source  response  set 


which  simplifies  t 


h2Z2  (Z|*  1)  -h,zl  (Z?*l) 


ratio- 

hJ>Jte-4T 

the  above  eguation  by  hlh2/hlh!  results  i 


which  is  written  i 


_(i-s)(sz2a(zaJi-i)-; 


MR  is  the  maximized  ratio  o 


e above  equation 


t (l-s)(sZ| (z|+l)  -Z?  (Z?»  1) ) 


RMR=maxa|-6s 
Setting  the  derivative  equal  to  2 


d(RMR)  . 


The  RMR  is  determined  from  s=0.577  (RMR=1.072).  The  point 
source  response  set  and  the  ray  used  to  calculate  the  RMR  are 
illustrated  in  Figure  4-8.  Since  the  RMR  was  searched  by  the 
ray  method,  the  solution  found  is  the  RMR  and  not  only  a local 
optimized  ratio. 

The  Utm^zation  of  the  RMR  for  the  Optimization  of  Detector 

As  mentioned  in  the  introduction,  the  RMR  can  aid  in 
assessing  the  performance  of  a detector  system.  Such  a 
performance  criterion  can  be  useful  for  both  the  operator  and 
designer  of  detector  systems.  While  the  calculational 
methodology  presented  in  this  chapter  is  sufficient  for 
assessing  the  RMR  of  a fixed  position  detector  system,  the 
more  complex  problem  of  the  design  a detector  system  and  the 
proper  detector  placement  of  a variable  position  detector 
system  is  not  addressed  in  depth  in  this  dissertation.  The 
technique  of  detector  system  design  should  include  the 
resulting  information  that  is  available  after  the  RMR  has  been 
calculated  in  order  to  attempt  and  improvement  in  the  RMR.  If 
the  line  intersection  method  is  used,  the  slopes  of  the  ray 
should  be  utilized  to  compare  the  relative  strengths  of  the 
detectors  that  contribute  to  the  RMR.  If  s is  formulated  as 
fftsf,,  then  s=f2/f2=tan0 . The  normalized  square  of  the 
magnitudes  of  the  vectors  f,  and  f,  that  contribute  to  the  ray 


defined  by  the  slope  s are 

norm (f?)  =cos20=cosJ(tan'1s)  and  norm(f.)  =sinJ(tan‘1s) . 
t}*ti 

The  normalized  square  of  the  magnitudes  give  the  relative 
strengths  of  the  detectors  that  contribute  to  the  RMR.  For 
example,  the  previous  example  gave  an  s of  0.577,  which 
results  in  norm (f,1) -0.75  and  normffj’l-O^S.  The  relative 
contributions  of  the  detectors  to  the  RMR  are  75%  for  f,  and 
25%  for  An  improved  measurement  configuration  could 
possibly  have  detector  #2  located  to  a position  that  is  less 
sensitive  to  positional  variations  of  the  source  positions. 

Very  close  to  the  source  container  the  spatial 
uncertainty  of  the  activity  contributes  significantly  to  the 
assessment  uncertainty.  As  the  measurement  location  increases 
in  distance  from  the  source  container,  the  statistical 
uncertainty  of  the  response  increase  resulting  in  the  optimal 
location  of  a detector  measurement  to  be  located  between  a 

infinite  distance  from  the  container.  For  simple 
detector/source  container  cases,  functional  expressions  can 
describe  the  relationship  between  the  RMR  and  the  location  of 
the  measurement.  Appendix  B contains  a development  of  an 
optimization  for  a one  and  a two  detector  configuration 
utilizing  a line  source  container. 


Figure  4-1.  Respective  supporting  planes  for  two  responses 
that  are  the  result  of  identical  distributions  and  different 
activities . 


Figure  4-2.  Arbitrary  supporting  planes 
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Figure  4-3 . Proper  supporting  planes  f 
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intersection  method 


Figure  4-5.  Line  intersection  method  for  complete  response 


Figure  4-6.  Method  of  searching  for  the  global  optimum 
utilizing  supporting  planes. 
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fixed  point  search. 
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CHAPTER 


ASYMPTOTIC  SOLUTIONS 


Increasing  number  of  detectors  to  infinity  does  not 
guarantee  that  the  assay  system  will  be  able  to  distinguish 
the  amount  of  activity  perfectly  (RMR=1) . As  the  number  of 
symmetric  detector  increases,  the  RMR  can  converge  to  some 
limiting  value  greater  than  one. 

Example  5-1.  Consider  line  source  container  located  on 
the  x axis  between  -0.5  and  0.5.  A number  of  symmetrically 
placed  detectors  are  located  on  the  midplane  of  the  source 
container  a distance  R from  the  origin  with  a response 


If  the  number  of  symmetric  detectors  is  increased  to  infinity, 
the  detector  system  cannot  distinguish  between  a unit  point 
source  located  at  the  origin  and  a point  source  of  strength 
(!?*+. 25)  /R2  located  at  x=.5.  If  the  source  container  of  the 
above  example  is  changed  to  a 55  gallon  drum,  increasing  the 
number  of  midplane  symmetric  detectors  would  aid  in 
distinguishing  source  positions  perpendicular  to  the  lateral 
axis  but  not  in  distinguishing  point  source  distributions 


along  the  axis.  Increasing  the  number  of  detectors  to  infinity 
in  this  case  would  cause  the  RMR  to  approach  the  asymptotic 
value  of  (R‘+.2S)/K‘. 

The  above  asymptotic  value  of  the  RMR  was  obtained  by 
inspection  of  the  source  container  and  the  detector  response 
equations.  The  asymptotic  value  can  be  approximated  by 
evaluating  the  RMR  as  the  number  of  detectors  becomes 
increasingly  large.  For  a generalized  detector  system,  a 
technique  for  the  direct  evaluation  of  the  asymptotic  RMR  has 
not  been  formulated  although  the  usefulness  of  such  a 
technique  has  been  mentioned  (Ben-Haim  1985b) . 

This  chapter  develops  a method  for  the  evaluation  of  the 
asymptotic  value  of  the  RMR  as  the  number  of  detectors 
increase  to  infinity.  The  analytical  solution  for  the 
asymptotic  RMR  is  solved  by  converting  the  set  in  detectors 
that  increase  to  infinity  to  an  integral.  The  following  method 
is  for  a set  of  symmetrical  detector  increasing  to  infinity 
along  with  a set  of  detectors  which  are  not  part  of  the 
symmetric  design.  Using  the  supporting  plane  method,  the  RMR 
for  the  N-l  symmetrically  placed  detectors  along  with  an 
additional  nonsymmetrically  placed  detector  is  written  as 


Since  the  first  N-l  detectors  are  symmetric,  the  corresponding 
weight  coefficients  are  equal.  Taking  into  account  that 


(5-3) 


(5-5) 


i,  j gndn*aKj  f„d 


As  demonstrated  in  the  following  example,  the  integral  can  be 
solved  analytically  by  converting  the  variable  n to  the 
spatial  parameters  of  the  problem. 

Example  5-2.  Consider  a source  confined  to  a circle  of 
radius  one  centered  on  the  origin.  The  N detectors  are  located 
symmetrically  on  a circle  of  radius  two.  Assume  that  the  Nth 
detector  is  fixed  at  (2,0).  Figure  (5-1)  illustrates  the  case 
of  four  symmetrically  placed  detectors  surrounding  a disk.  The 
source  location  that  corresponds  to  the  minimum  detector 
response  for  any  number  of  detectors  greater  than  one  is 
located  at  the  center  of  the  disk  (0,0).  The  maximized 
detector  response  for  the  nth  detector  located  at 


(5-7) 


The  asymptotic  limit 
parameter  is 


The  solution 


( i dfl 

N J 5-4COS0 

2s  j 1 dB 


(5-8) 


calculation 


to  the  above  equation  is  4/3.  A numeric 
eighteen  symmetrically  placed  detectors  results 


The  above  solution  considers  the  multiple  responses  as 
multiple  detectors.  Other  perspectives  can  be  taken  which  will 
result  in  the  same  formulation  of  the  problem.  For  time 
invariant  problems,  a single  detector  at  multiple  positions 
can  be  considered  identical  to  multiple  detectors  at  fixed 
locations.  Another  alternative  would  be  to  have  a single 
detector  remain  stationary  and  have  the  source  container  move 
to  another  position  or  rotate  about  an  axis.  Such  a rotation 
has  been  previously  considered  (Martens  and  Fliss  1990)  as  a 
means  of  reducing  the  uncertainty  assessment. 

As  the  number  of  the  symmetric  detectors  approach 
infinity  the  resulting  "infinite-symmetric"  detector  is 
represented  by  an  integral  which  is  equivalent  to  a single 
detector  with  the  same  geometry  that  the  infinite  number  of 
point  detectors  form.  In  the  above  example,  the  infinite 
number  of  detectors  surrounding  the  disk  is  equivalent  to  a 
ring  detector.  In  the  case  of  the  above  example,  the  problem 
can  be  redefined  as  a rotating  source  container  and  a detector 
located  somewhere  on  the  circle  r=2 . The  response  function  for 


circular 


container  rotating 


with  a point  source  located  at  ( r, 0)  is 

f, - f —= — -c(8. 

{ R1+z1-RzcosQ 

Multiplying  and  dividing  by  Jf3  gives 


(5-9) 


(5-10) 


A table  of  definite  integrals  (Spiegel  1968,  15.47)  provides 
the  following  solution 


The  minimum  detector  response  occurs 
response  occurs  at  the  edge  of  the  sou 
RMR  of 


(5-11) 

r=0  and  the  maximum 
container  giving  an 


(5-12) 


which  for  a circular  source  container  of  radius  one  and  the 
detector  placed  at  R= 2 results  in  an  RMR  of  4/3. 

The  asymptotic  limit  of  the  RMR  can  assist  the  assessor 
in  developing  a measurement  strategy  that  maximizes  the 
efficiency  of  the  detector  measurements  while  minimizing  the 
time  and  effort  necessary  to  achieve  those  measurement.  The 
RMR  for  an  infinite  number  of  detector  in  the  above  example 
was  4/3.  This  asymptotic  limit  may  be 


achieved  by  a rotation 


of  the  container  or  a ring  detector.  Four  symmetrical 
measurements  around  the  disk  achieve  an  RMR  of  1.51  and  eight 
symmetric  detector  achieve  an  RMR  of  1.35.  Therefore,  the 
assessor  can  determine  whether  the  time  and  expense  of 
increasing  the  number  of  detector  measurements  is  worthwhile. 


symmetric  detectors  surrounding 


CHAPTER 


The  application  of  Euclidean  space  to  the  problem  of 
assessment  of  radioactive  material  in  the  presence  of  spatial 
uncertainty  has  achieved  limited  use  as  both  a research  and  an 
industry  tool.  More  common  is  the  use  of  methods  such  as 
segmentation  or  rotation  of  the  container  and  a linear 
combination  of  responses.  These  methods,  while  able  to  reduce 
the  uncertainty,  do  not  have  the  flexibility  to  adjust  to 
modifications  in  container  geometry,  detector  positions  or 
number  of  detectors.  As  a corollary,  the  above  methods  do  not 
allow  for  an  expansion  to  concepts  of  optimization  and 
feedback  control. 

Without  application  of  Euclidean  space,  the  equations  for 
interpretation  and  optimization  rely  on  factors  such  as  the 
symmetry  the  detector  source  configuration,  assumed  spatial 
distributions  of  the  source  and  a loss  of  information  (through 
averaging  or  a linear  combination  of  the  multiple  responses) 
to  allow  for  a solution  to  the  problem.  The  application  of 
dimensional  space  allows  for  a mathematical  framework 
independent  of  container  geometry,  measurement  position  and 
source  distribution.  This  dissertation  has  expanded  the 


application  of  Euclidean  space  to  develop  solutions  for 
problems  that  have  previously  been  unsolved,  such  as  an 
analytical  solution  to  a asymptotic  limit,  the  utilization  of 
nonconvex  response  sets,  the  global  optimization  of  the  RMR 
and  the  interpretation  of  multiple  measurements.  The  results 
are  taken  a step  further  by  combining  the  interpretation 
equation  with  the  RMR  equation  to  develop  a feedback 
optimization  procedure  where  the  results  of  the  previous 
measurements  aid  in  determining  the  position  of  a future 

A possible  reason  for  the  limited  use  of  a vector  space 
approach  for  multiple  measurements  is  the  lack  of  tools 
necessary  for  the  assessment  of  the  radioactive  material 
measured.  The  formulations  developed  in  this  dissertation, 
especially  the  interpretations 1 method,  demonstrate  that  the 
application  of  a multidimensional  space  provides  generalized 
solutions  to  a variety  of  assessment  problems. 

Since  a discrepancy  exists  between  the  potential  of  the 
application  of  vector  space  and  the  number  of  references 
available  on  the  subject,  there  are  other  significant 
accomplishments  incorporated  in  the  dissertation.  The 
inclusion  of  simplified  examples,  the  "filling  in"  of 
information  assumed  in  other  literature,  the  "walking  through" 
of  results,  the  emphasis  on  unifying  concepts  and 
interpreting  problems  from  multiple  perspectives  should  allow 


this  dissertation  to  serve  as  a reference  and  tutorial 
would-be  practitioners  and  researchers. 


4-S--4 


f=(ag*ll-a) h)  where 


(A-8) 


A set  F is  convex  if  for  each  pair  of  points  (g,,g3t  - - - , gN)  and 
(h„h3, . . • ,hN)  in  the  set  C,  the  line  segment  connecting  these 
points  is  also  in  the  set  C.  A convex  combination  of  M points, 
g,h,  is  defined  as  the  set  of  points  (c„c„  such 


c=ayg'a2h’ . . . •aHk  where  £ a„=l  and  all  a„21 . (A-9) 

If  the  number  of  points  M approaches  infinity,  the  constraint 
of  the  coefficient  is 


faadn= 1 and  all 

convex  combination  of  points  of  a : 
is  represented  as 

C-{c:c=j\t„fdn  such  that  j" a, 


(A— 10) 

let  C which  are  the 
composed  of  points 

and  «„*!) . (A-ll) 


The  relationship  between  a convex  set  and  a convex  combination 
of  points  is  as  follows:  a set  is  convex  if  and  only  if  every 
convex  combination  of  points  in  S lies  in  S (Ben-Haim  1985a) . 

A convex  hull  of  a set  s is  defined  as  the  intersection 
of  all  the  convex  sets  which  contain  S and  is  denoted  as  the 
convex  hull  of  S.  A convex  hull  is  the  smallest  convex  set 
that  contains  S.  In  addition,  S is  convex  if  and  only  if  S 
equals  the  convex  hull  of  S (Lay  1982) . The  convex  hull  of  S 


consists  of  exactly  all  convex  combinations  of  elements  in  the 
set  S.  A supporting  plane  of  a set  S is  defined  as  the  plane 
that  has  at  least  one  point  in  common  with  the  set  s and  is  a 
maximum  or  minimum  distance  from  the  origin.  The  upper 
supporting  plane  is  a maximum  distance  from  the  origin  and  the 
lower  supporting  plane  is  a minimum  distance  from  the  origin. 
A supporting  plane  of  the  set  S is  written  as 

maxfor  min)MS<a,h>.  (A-12) 

A supporting  plane  in  2 dimension  space  is  a line.  Figure  A-3 
illustrates  an  example  of  an 
for  a two  dimensional  convex 


upper  and  lower  supporting  plane 


A 


ft 


Directional 


i 


intersected 


from  the  origin. 


supporting  planes  for 


APPENDIX  B 

ADAPTIVE  ASSAY  WITH  STATISTICAL  AND  SPATIAL  UNCERTAINTY 


so  far  been  ignored,  although  the  formulations  for  the  RMR 
that  include  statistical  uncertainty  have  been  derived  (Ben- 


detector  responses.  Without  the  incorporation  of  statistical 
uncertainty  the  ideal  detector  configuration  will  frequently 


quantities 


(A,(m*5)(:u)  > (h.m tv). 


(B-2) 


If  o is  the  statistical  standard  deviation,  then  the  RMR  is 
based  on  the  equation 

(h,  (m»8)  tu)-bo((h,  (m-8)  tu))  >(h,mtv)*bo((b,  me  v»  (B-3) 
where  b is  a positive  number  indicating  the  level  of 
statistical  confidence  required.  Solving  the  above  equation 
for  the  RMR  results  in 

RMR  ■ =(.JgMR*V1+T*v)2  (B-4) 


T-  ba«h,  v))  ^ ba«h,  v)) 


(B-5) 

Note  that  the  amount  of  source  material  present,  which  is 
to  be  interpreted  from  the  detector  measurements,  is  required 
for  the  determination  of  the  RMR.  If  the  amount  of  mass  can  be 
approximated,  the  detectors  can  be  optimized  prior  to 
measuring  the  sample.  Since  the  amount  of  radioactive  material 
is  the  variable  that  is  being  assessed,  however,  the  mass  may 
not  be  known  with  any  degree  of  accuracy.  The  assessment  of 
the  amount  of  radioactivity  present  from  an  initial  detector 
configuration  provides  information  in  order  that  the  locations 
of  the  detectors  can  be  changed.  The  new  locations  of  the 
detectors  will  result  in  a lower  RMR  for  the  detector  system 
and  possibly  a reduced  range  of  the  interpreted  amount  of 
radioactive  material  present.  The  initial  m 


incorporated  into  the  readjusted  measurements  by  forming  an 
N+R  detector  space  where  N is  the  number  of  initial 
measurements  and  R is  the  number  of  relocated  detector 

A brief  remark  was  previously  presented  (Ben-Haim  1985b, 
pg.  305)  where  the  design  method  of  minimizing  the  RMR  is 
performed  by  adjusting  the  Nth  detector  position  based  on  data 
from  the  previous  N-l  detector  positions.  This  suggestion 
utilizes  the  previous  N-l  detector  measurements  taken, 
incorporates  these  measurement  within  an  N detector  space  and 
attempts  an  optimization  of  the  Nth  detector  response.  For 
example,  a single  detector  interprets  the  amount  of 
radioactive  material  from  the  detector  response.  The 
assessment  of  mass  indicates  that  the  detector  position  should 
be  changed.  This  new  detector  location  provides  an  RMR  that  is 
lower  than  the  RMR  from  the  previous  detector  location. 
However,  if  the  source  configuration  is  time  independent,  the 
previous  detector  response  and  the  optimized  detector  response 
can  be  considered  as  a two  detector  set-up.  The  RMR  for  this 
two  detector  set-up  will  be  lower  or  equal  to  the  optimized 
single  detector  set-up. 

The  following  example  problem  confines  the  radioactivity 
to  a line  of  unit  length  along  the  x-axis.  The  optimal 
location  of  the  detector  is  found  by  utilizing  the  information 
from  the  first  measurement  and  creating  a two  dimensional 
response  space.  The  detector  response  located  a distance  R,  to 


*—0.5;  therefore,  the  FMR  is 

—tim 


The  average  of  the  two  responses  is  2.223.  Using  the  above 
equation,  a single  detector  should  be  located  at  £,=1.721. 

In  order  to  incorporate  the  previous  measurement  of  0.889 
of  a future  measurement,  a two  detector  system  is  established 
with  the  first  detector  fixed  at  x=1.5.  The  optimum  location 
of  the  second  detector  may  not  be  at  £,=1.721.  Assume  that  the 
readjusted  detector  is  located  somewhere  on  the  positive  pr- 
axis. The  detector  response  is 


The  function  will  be  attempted  to  be  classified  as  concave  or 
convex  by  solving  for  the  second  derivative.  The  response  of 
detector  1 in  terms  of  detector  2 is 

f,  = 1 . 

h'^-i-s))  (b_i2) 

The  derivative  is 


dfi  (/^♦1.5)-l),-2f^/II<£,tl.S)  -Ilj-^IJVI.S)] 

"57‘"  (/?!(£,  .1.5)-!)' 

Simplifying  the  expression  gives 

<**2  _ • 5)-l)M^ • 5)  -1)(^ (£,+1  ■ 5) ) (B-14J 

dt'  l.S)-!)1 


df;,  -3(/r^i.5)-i) 
df-  (/?!(«, ‘1.5J-1)1 


The  numerator  of  the  second  derivative  is 


(B-15) 


md2f,  / 3 (Rj-t-l . 5) 

dfi  "l  2V^  . 


Factoring  out  gives 


— g ■(  ’ )•«■»  «-»■> 

which  is  greater  than  zero;  therefore,  the  function  is  convex 
and  the  upper  supporting  plane  will  intersect  the  endpoints 
and  the  lower  supporting  plane  will  be  tangent  to  the  point 
source  response  function. 

The  endpoints  of  the  point  source  response  function  occur 
at  x=-.5  and  x=.S.  The  slope  of  the  line  connecting  these 


f,(.5)-f,(-.5)  ,a,  -R, 

The  derivative  of  f2  with  respect  to  f,  in  terms  of  x is 

d^-df,dx_  2 (1.5+x)3_  (1.5+x)3 

df,  dx  df,  (Rj-x)1  1 (Rj-x)3 


(B-20) 


tangent  at  x=x,  where 


the  equation 


_/  8\ 

(Kj-X.)1  I 3/  (l?2~ . 5) 2 (i?2+5)  2 * 


Solving  for  x,  gives 


(B-22) 


The  RMR  with  statistical  uncertainty  is 


(B-24) 


3(l^-.5)>(i^+.5)a(l.S 


*U-x.)  (B.26) 

16  Ra^5 + 2—^— 

3^R1-.S)1(R1*.S)1(1.S*X.)^ 

The  above  equation  is  dependent  only  on  R;  and  m.  Using  one 
detector,  the  amount  of  activity  present  was  calculated  to  be 
between  0.889  and  3.556,  an  average  of  2.223.  The  optimal 
location  for  the  detector  is  found  by  substituting  an 
estimated  value  for  m,  setting  the  derivative  of  the  RMR  with 
respect  to  R,  equal  to  aero  and  solving  for  R,. 

Since  the  previous  detector  response  of  f,-0. 889  will  be 
used  for  the  interpretation  of  the  amount  of  radioactivity 
present,  the  above  solution  determined  the  RMR  without  regard 
to  the  information  from  the  previous  detector  response.  The 
response  of  detector  / 1 limits  the  assessment  uncertainty 
present.  If  the  ray  method  is  utilised  for  determining  the  RMR 
of  the  above  example,  the  ray  intersecting  the  response  set  is 
limited  to  (0.889, f.)  where  f2  is  confined  to  the  complete 
response  set  multiplied  by  the  estimated  amount  of 
radioactivity.  The  restriction,  f,=0.889,  places  the  point  of 
intersection  of  the  ray  at  (0.889,f;)  where  f,  could  result  in 
the  point  being  on  the  boundary  of  the  complete  response  set 
or  a point  within  the  complete  response  set.  The  equation  of 
the  point  source  response  of  detector  f 2 is  dependent  only  on 
the  location  of  the  detector  and  the  source  amount  present 


fa. £i£ 


(B-27) 


For  purposes  of  calculating  the  RMR  the  value  of  m can  be 
considered  one  since  the  RMR  of  the  identical  regardless  of 
the  amount  of  radioactive  material  assumed.  Assuming  f, =0.869 
and  m=l,  the  equation  of  the  ray  that  intersects  the  point  is 


The  equation  of  the  line  that  intersects  the  two  endpoints  of 
the  point  source  response  set  is 


The  intersection  of  this  line,  which  bounds  the  complete 
response  set,  with  the  ray  is 


The  intersection  with  the  point  source  response  set  occurs  at 


(B-29) 


■943*j*.414  l.S>-l)*. 


(B-31) 


*1  (*,*1 . 5)  s-2  (*j*l . 5)  - . 943R;,* . ! 

quadratic  equation  results  in 


(B-32) 


/,=2,/I7(*,< 


^«/,  (*,♦!■  5)^-4  (*,♦!■ 5)  J(-.943*,+  .l 
2(R,*1.5)J 


d simplifying  gives 


l±/.943*,*-‘ 

(*,♦1.5) 


The  above  analysis  is  incomplete  in  that  it  assumes  that 
the  point  (0.889,/,)  is  located  on  the  point  source  response 
set.  A completed  analysis  will  include  that  (0.889,/,)  is  part 
of  the  other  boundary  of  the  complete  response  set.  The  RMR 
for  the  problem  is  between  the  two  calculated  RMRs  from  the 
boundary  points.  The  optimized  detector  location,  in  addition, 
will  have  a range  of  calculated  values. 

The  concept  of  the  time  variable  t has  been  ignored  (Ben- 
Haim  1985b) . If  total  measurement  time  available  is  T,  and  the 
first  measurement  used  t,  units  of  time,  then  T-t,  units  are 
available  for  the  next  measurement.  The  uncertainty  resulting 
from  a second  measurement  for  the  entire  time  remaining,  T-t, 
may  be  too  large.  The  determination  the  number  of  detector 


the  partitic 


tine 


those  measurements  in  order  to  achieve  a desired  level  of 
certainty  remains  a reletively  unexplored  topic. 
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